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I.  INTRODUCTION 


Projectile  impact  shock  initiation  of  high  explosives  has  long  been  a 
subject  of  interest  in  the  energetic  materials  community.  Considerable 

experimental  data  has  been  generated  over  the  years.  Numerical  modeling  of 
projectile  impact  shock  initiation  for  comparison  with  experiments  has  been 

4 

reported  in  at  least  one  case.  However,  only  the  predicted  critical  velocities 
and  no  detailed  analysis  of  the  flow  fields  revealed  by  the  computations  were 

5 

presented.  In  another  numerical  study,  Mader  and  Pimbley  modeled  the 
initiation  of  explosives  due  to  the  impact  of  shaped  charge  jets  using  the  same 
computer  code  used  in  the  present  report. 

This  report  concerns  our  numerical  modeling  of  the  projectile  impact  shock 
initiation  of  composition-B  (oomp-B).  We  have  considered  both  bare  and  covered 
charges  impacted  by  cylindrical  steel  projectiles  using  the  Los  Alamos  2DE 
code.  We  have  examined  the  flow  fields  in  some  detail  and  compared  predicted 
critical  velocities  with  published  experimental  values.  Our  earlier  work  on 

this  subject  was  found  to  have  made  use  of  insufficient  artificial  viscosity.® 
This  has  been  corrected  in  the  present  report  and  the  work  extended  to  a  number 
of  different  areas. 


lD.  C.  Slade  and  J.  Dewey,  " High-Order  Initiation  of  Two  Military  Explosives 
by  Projectile  Impact,  ”  Ballistic  Research  Laboratory  Report  No.  1021 ,  July 
1957.  AD  145868 

2S.  M.  Brown  and  E.  G.  Whitbread,  "The  Initiation  of  Detonation  by  Shock 
Waves  of  Known  Duration  and  Intensity,  "  Les  Ondes  de  Detonation ,  C.N.R.S. 

No.  109,  pp.  69-80,  Paris,  1962.  ■ 

3i.  A .  Roslund ,  J.  W.  Watt,  and  N.  L.  Colebum,  " Initiation  of  Warhead 
Explosives  by  the  Impact  of  Controlled  Fragments,  I  Normal  Impact,  "  Naval 
Ordnance  Laboratory  Technical  Report  NOLTR-73-124,  August  1974. . 
kK.  '.  Bahl,  R.  C.  Vantine,  and  R.  L.  Weingarts,  " The  Shock  Initiation  of  Bare 
and  Covered  Explosives  by  Projectile  Impact,  "  Seventh  Symposium  ( International ) 

on  Detonation,  June  1981,  pp.  325-335.  '  ■  . 

5 c.  l.  Mader  and  G.  R.  Pimbley ,  "Jet  Initiation  and  Penetration  of  Explosives, 
Journal  of  Energetic  Materials,  Vol.  1,  No.  1,  1983,  pp  1*44.  _  ^  ^ 

tf.  X.  Huang',  J.  Starkenherg'," and  A.  L.  Arhuckle,  "A  Numerical  Study  of  Shock  ^ 
Initiation  of  Composition-B  by  High  Speed  Impact  of  Small  Steel  Projectiles, 
BRL  Report  to  be  published. 


II.  BRIEF  DESCRIPTION  OF  20E 

The  2DE  code'  is  s  two-dimensional,  reactive  hydrodynamic  computer  code 
which  makes  use  of  the  equations  of  motion  in  Eulerlan  form.  It  incorporates 
the  HOM  equation  of  state  andt  most  Important  for  our  application,  the  Forest 
Fire  model  for  shock  initiation  of  high  explosives.  In  our  computations  we  used 
an  elastic-plastic  constitutive  model  to  account  for  the  behayior  of  steel. 


III.  PROJECTILE  IMPACT  SHOCK  INITIATION  OF  BARE  COHPOSITIOM-B 
Geometry  and  Computational  Considerations 

The  axisymmetric  geometry  used  in  the  bare  charge  projectile  impact 
computations  is  illustrated  in  Figure  1.  Ve  have  considered  cylindrical  steel 

projectiles  of  unit  aspect  ratio  U/d)  since  Brown  and  Whitbread  have 
demonstrated  that  different  aspect  ratios  do  not  produce  different  critical 
velocltes  for  shock  initiation,  except  in  the  ease  of  thin  discs  (£/d<1/4). 
Computations  were  made  for  projectile  diameters,  d,  of  4,  5,  8,  10,  12,  15  and 
18  mm.  Sufficient  target  material  is  provided  when  the  length  and  diameter  of 
the  explosive  charge  are  each  three  times  the  corresponding  projectile 
dimensions. 

Ve  set  up  these  impact  problems  for  20E  calculation  with  axlsymmetrlc  grids 
as  summarized  in  Table  1.  Here  Ar  is  the  radial  cell  dimension,  I  the  number 
of  cells  slong  the  radial  axis,  az  the  axial  cell  dimension,  J  the  number  of 
cells  along  the  axis  of  symmetry.  At  the  time  step  for  each  computational  cycle, 
and  N  the  total  number  of  cycles  to  be  completed. 


TABLE  l 

Input  Oats  for  2DE  Computational  Grids  -  Bare  Charges 


d 

(mm) 

Ar,A* 

(mm) 

I 

J 

At 

(ms) 

N 

4 

0.200 

40 

90 

0.005 

400 

5 

0.250  _ 

_ 40 

70 

0.006 

400 

8 

0.207 

60 

105 

0.006 

500 

10 

0.340 

60 

105 

0.010 

500 

12 

0.400 

60 

105 

0.008 

500 

15 

0,500 

60 

105 

0.010 

50U 

18 

0.600 

60 

105 

0.010 

400 

7J.D.  Kershter  and  C.  L.  Made?,  "2DE,  A  Tuo-Dimensional  Continuous  Eulerlan 
Hydrodynamic  Code  for  Computing  Multicomponent  Reactive.  Hydrodynamic  Prob¬ 
lems,  "  Los  Alamos  Scientific  laboratory  Report  LA-4848,  March  1972 . 
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Axisymmetric  Geometry  Used  in  Projectile  Impact  Computations. 
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Results 


Flow  Field  Observations.  A  number  of  graphical  representations  of  our 
numerical  results  are  available.  The  sequence  of  events  in  projectile  impact 
shock  initiation  is  most  clearly  illustrated  in  the  series  of  mass  fraction 
contour  plots  of  Figure  2.  The  mass  fraction  varies  from  one  to  it.ro  as 
chemical  reaction  in  the  explosive  runs  to  completion.  These  plots  show  results 
for  the  impact  of  a  10  mm  diameter  projectile  at  1.1  ka/s,  just  above  the 
critical  velocity.  The  position  of  the  impact  shock  is  also  shown.  At  1.0  ms 
after  impact,  burning  is  observed  throughout  the  region  between  the  shock  and 
the  projectile,  but  detonation  has  not  yet  begun.  Detonation,  which  may  be 
recognized  by  the  close  spacing  of  the  contour  lines,  is  first  observed  to  break 
out  after  the  shock  has  propagated  some  distance  from  the  Impact  point.  The 
detonation  then  spreads  outward  along  the  shock  and  is  well  established  by 
1.5  ms.  It  also  propagates  back  toward  the  projectile  through  the  partially 
reacted  material. 

Determination  of  Critical  Impact  Velocity  and  Comparison  with  Experimental  Data. 
By  varying  impact  velocity  in  computations  of  this  type  we  were  able  to 
determine  the  critical  velocity  as  a  function  of  projectile  diameter.  In  Figure 
3  we  have  plotted  the  2DE  results  together  with  the  experimental  results  of 

*  1  3 

Dewey  and  Slade  as  well  as  the  Jacobs-Roslund  empirical  formula  .  The 

agreement  is  excellent,  with  the  2DE  go  and  no-go  predictions  all  bracketing  the 

empirical  curve.  We  have  also  included  the  data  for  ISL  comp-B  (65/35)  reported 

by  Moulard  at  the  last  detonation  symposium.**  .This  explosive  is  reported  to  be 
generally  less  sensitive  than  U.S.  comp-B  and  particularly  less  sensitive  when 
impacted  by  small  projectiles! 

We  observed  a  different  mode  of  critical  shock  initiation  at  the  smaller 
diameters.  When  the  projectile  diameter  is  small  and  the  impact  velocity  is 
high^ detonation  appears  almost  immediately  on  impact  as  shown  in  Figure  4  at 
0.3  ms  after  the  impact  of  a  4  am  diameter  projectile  at  1.7  km/s.  In  this  case 
detonation  continues  to  propagate  and  a  considerable  amount  of  explosive  has 
been  consumed  by  1.2  ms.  On  the  other  hand,  when  the  impact  velocity  is 
reduced  to  1.6  km/s,  the  detonation  that  breaks  out.  immediately  does  not 
continue  to  propagate  but  is  quenched  by  the  action  of  following  rarefactions  as 
shown  in  Figure  5.  Thus,  the  mass  fraction  contour  lines  begin  to  spread  out  by 
0.7  ms.  Little  or  no  progress  is  made  between  1.0  ms  and  1.5  ms  as  the 
detonation  dies  out  leaving  a  bubble  of  detonation  products  in  its  wake.  It 
should  be  noted  that  these  detonations  are  not  overdriven  as  were  those  computed 

by  Nader  and  Pimbley^  for  shaped  charge  jet  impact.  An  overdriven  detonation  in 
comp-B  produced  by  the  impact  of  a  steel  projectile  would  require  an  impact 
velocity  exceeding  approximately  2.8  km/s.  Actually,  full  CJ  pressure  is  never 
achieved  in  these  detonations.  Figure  6  shows  a  series  of  pressure  and  mass 


8ff.  Moulard ,  "Critical  Conditions  for  Shook  Initiation  of  Detonation  by  Smalt 
Projectile  Impact, "  Seventh  Symposium  (International)  on  Detonation,  June 
1981,  pp.  816*824. 
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fraction  profiles  along  the  axis  of  symmetry  at  various  tiaes  for  the  4  aa 
projectile  impact  at  1.6  ka/s.  While  the  aass  fraction  drops  rapidly  to  zero 
the  pressure  never  rises  above  about  23  GPa. 

In  order  to  assess  the  relationship  between  the  2DE  predictions  of  critical 

2 

velocity  and  the  critical  p  t  concept,  we  aade  a  series  of  computations  in  which 
the  Forest  Fire  model  was  deactivated  and  the  explosive  treated  as  an  Inert 
material.  By  observing  the  pressure  history  of  the  target  explosive  adjacent  to 

the  lapact  point  we  were  able  to  calculate  dt.  Computations  corresponding  to 

our  highest  subcritlcal  impact  velocities  for  projectile  diameters  ranging  from 
5  am  to  18  mm  were  made.  The  results,  suansrlzed  in  Table  2,  show  that,  while 

peak  shock  pressure  decreases  with  decreasing  lapact  velocity,  ^p2dt  remains 

fairly  constant  along  the  Initiation  threshold.  Thus,  the  critical  p2t  concept 

appears  consistent  with  Forest  Fire. 

Table  2,  Response  of  "Inert”  Composition-B  to  Critical  Frojectile  Impact 


Projectile 

Diameter 

(am) 

Impact 

Velocity 

(ka/s) 

Peak  Shock 
Pressure 

(GPa) 

/p2dt 
(GPa2-  ms) 

5 

1.4 

9.6 

32 

8 

1.1 

7.0 

32 

10 

1.0 

6.0 

31 

12 

0.9 

5.3 

28 

15 

0.8 

4.5 

28 

18 

0.7 

3.8 

25 

Shock  to  Detonation  Transition  Paths.  The  Forest  Fire  Model  is  based  in  part  on 
the  Single  curve  buildup  hypothesis.  Thus,  the  Pop-plot  is  interpreted  as 
describing  the  process  of  buildup  to  detonation  in  the  shock  pressure  -  distance 
to  detonation  plane.  This  is  true  at  least  for.  the  planar  geometries  in  which 
single  curve  buildup  has  been  observed.  In  an  earlier  numerical  study  using 
2DE,  we  had  occasion  to  consider  the  planar  problem  .art sing  when  a  flyer  plate 

e 

of  sufficient  lateral  extent  strikes  a  eomp-B  target.  Once  the  distance  along 
the  axis  of  symmetry  at  which  steady  state  detonation  first  appears  has  been 
determined,  the  progress  of  shock  buildup  toward  detonation  as  a  function  of 
distance  of  run  to  detonation  can  be  compared  to  the  Pop- plot.  This  has  been 
dona  for  the  problems  of  10  mm  thick  flyer  plates  striking  comp-B  targets  at  0.6 
and  0.7  ka/s  in  Figure  7.  The  results  indicate  that  2DE  reproduces  the  single 
curve  buildup  phenomenon  in  the  planar  geometry.  The , pro jectile  impact  data, 
however,  does  not  appear  to  produce  a  single  curve  buildup  along  the  axis  of 
symmetry  as  illustrated  in  Figure  8. 


V.  Starkenberg ,  7.  K.  Huang  t  and  A .  L.  Arbuokle ,  *4  Tuo~Dimensional  Numerical 
Study  07  Detonation  Propagation  Between  Munitions  by  Means  of  Shook  Initiation , " 
BRL  Report  ARBRL-TR-02522,  September  198$.  ADA  T  33680 
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Figure  2 .  Sequence  of  Mass  Fraction  Contour  Plots  for  the  Supercritical 
Impact  of  a  10  mm  Diameter  Steel  Projectile  at  1.1  ka/s  against 
Bare  Composition-B 
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Figure  2*  (continued) 


13 


i 


0 


Figure  3.  Critical  Impact  Velocity  a*  a  function  of  Projectile  Diameter  -  Comparison 
2DE  Predictions  with  Published  Experimental  Data  for  Bare  Composition-B. 
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Figure  4.  Sequence  of  Mass  Fraction  Contour  Plots  for  the  Supercritical 
Impact  of  a  4  m  Diameter  Steel  Projectile  at  1.7  km/s  against 
Bare  Coaposition-B 
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Figure  5.  Sequence  of  Mass  Fraction  Contour  Plots  for  the  Subcritical 

Impact  of  a  4  mm  Diameter  Steel  projectile  at  1.6  km/s  against 
Bare  Composition-B 
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Figure  5.  (continued) 
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Figure  6.  (continued) 
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Impact  against  Bare  Composition- 


Special  Geometric*.  At  the  Seventh  Symposium  on  Detonation,  Houlard  presented 
som  interesting  experimental  results  in  which  ISL  eoep-B  was  impacted  by 

8 

projectiles  of  circular,  annular  and  rectangular  cross-section.  Although  the 
projectiles  were  designed  to  yield  the  same  overall  shock  duration,  he  observed 
considerable  differences  in  the  critical  velocities  produced  by  each. 
Specifically  he  found  that  the  cylindrical  projectile  required  the  highest 
critical  velocity  (about  2  km/s),  the  atmular  projectile  required  the  lowest 
critical  velocity  (less  than  1  km/s),  and  the  rectangular  cross-section  produced 
an  intermediate  critical  velocity.  He  sought  to  explain  this  by  Introducing  a 
critical  area  concept.  We  were  interested  in  determining  whether  the  classical 
shock  initiation  concepts  incorporated  in  2DE  could  explain  these  observations 
without  recourse  to  additional  concepts.  The  circular  end  annular 
cross-sections  could  be  represented  exactly.  Indeed,  the  circular  cross-section 
computation  had  already  been  completed.  The  impact  of  a  rectangular 
cross-section  projectile  is,  strictly,  a  three  dimensional  problem,  but  we 
represented  it  by  the  impact  of  a  slab  of  infinite  breadth  having  the  thickness 
of  the  smallest  dimension  (5  mm)  of  the  projectile  used  in  the  experiments.  The 
experimentally  and  numerically  determined  critical  velocities  are  listed  in 
Table  3.  The  results  show  that  the  classical  concepts,  on  which  the  Forest 
Fire  Model  is  built  are  sufficient  for  a  qualitative  explanation  of  the  Moulard 
observations.  The  principal  reason  for  the  absence  of  quantitative  agreement  is 
the  different  formulation  of  ISL  comp-B  and  its  reported  lower  senstivlty  to 
small  diameter  projectile  impact. 

Table  3.  Comparison  of  Moulard  Experiments  with  2DE  Simulation. 

Critical  Velocity  (km/s) 

Projectile 


Cross-section 

Moulard  Experiments 

2DE  Predictions 

(65/35  comp-B) 

(60/40  comp-B) 

5  mm  @ 

1.95  -  2.02 

1.40  -  1.50 

5  mm  i.d. 

IS  mm  o.d. 

1.06  -  1.11 

0.80  -  0.90 

mm  x  11  mm 

1.33  -  1.42 

1.10  -  1.15 

Some  of  the  reasons  for  the  variations  in  critical  velocity  Mi  projectile 
cross-section  become  apparent  when  we  observe  the  flow  fields  pro  ^ced.  Figure 
9  shows  a  series  of  mass  fraction  contour  plots  for  the  impact  of  the  annular 
projectile.  The  relatively  low  impact  velocity  produces  no  immediate  reaction 
adjacent  to  the  Impact  point.  However,  at  a  later  time,  shock  reflection 
(probably  Mach  stem  formation)  at  the  axis  of  symmetry  produces  higher  pressures 
than  the  circular  cross-section  projectile  impact  at  the  same  velocity.  Thus, 
the  detonation  is  observed  to  break  out  along  the  axis  of  symmetry  at  lower 
impact  velocities.  Of  perhaps  greater  Interest,  then,  is  the  difference  between 


the  circular  and  rectangular  cross-section  results.  The  only  major  dii Terence 
here  is  the  rate  at  which  the  rarefaction  quenches  the  incipient  detonation. 
Remember  that  in  the  small  diameter  case  we  observed  immediate  detonation  which 
was  then  ouenched  by  the  action  of  the  following  rarefaction.  The  results  for 
the  rectangular  cross-section  are  shown  in  Figures  10  and  11.  Figure  10  shows 
mass  fraction  contours  for  a  supercritical  impact  at  1.25  km/s.  Detonation 
arises  as  a  result  of  shock  to  detonation  transition.  In  the  subcritical  impact 
at  1.1  ka/s  in  Figure  11,  no  detonation  occurs.  It  remains  to  be  determined 
whether  this  strong  effect  on  critical  velocity  is  manifested  with  projectiles 
of  larger  dimensions  for  which  simple  shock  to  detonation  transition  occurs  for 
the  cylindrical  projectiles  also. 

IV.  PROJECTILE  IMPACT  SHOCK  INITIATION  OF  COVERED  COMPOSITION- 2 
Geometry  and  Computational  Considerations 

We  have  also  addressed  the  related  problea  of  projectile  impact  of  covered 
comp-B  by  introducing  a  steel  plate  of  thickness,  h,  between  the  projectile  and 
the  explosive.  Projectile  diameters  of  6,  8,  and  10  am  and  cover  plates  of  1/8, 
1/4,  1/2,  and  3/4  as  thick  as  the  diameter  in  each  case  were  considered.  Other 
geometrical  considerations  are  as  described  for  the  bare  charge  problem.  The 
axisymaetric  computational  grids  are  described  in  Table  4. 

Table  4.  Data  for  2DE  Computational  Grids  -  Covered  Charges 


d 

(mm) 

h/d 

Ar.AZ 

(mm) 

I 

J 

At 

(ms) 

N 

6 

1/8 

0.15 

45 

100 

0.004 

400 

6 

1/4 

0.15 

45 

100 

0.004 

600 

6 

1/2 

0.15 

45 

100 

0.004 

650 

6 

3/4 

0.15 

45 

100 

0.004 

700 

8 

1/8 

0.20 

45 

100 

0.005 

450 

8 

1/4 

0.20 

45 

100 

0.005 

600 

8 

1/2 

0 .20 

45 

110 

0.005 

8W 

8 

3/4 

0.20 

45 

,120 

0.005 

800 

10 

1/8 

0.25 

45 

100 

0.005 

500 

10 

1/4 

0.34 

55 

125 

0.008 

700 

10 

1/2 

0.20 

55 

140 

0.005 

750 

10 

_ 2/1— 

0.20 

55 

140 

0.005 

800 

Results 

Flow  Field  Observations.  Typical  results  are  shown  in  the  mass  fractijn  contour 
plot  sequence  of  Figure  12.  This  is  for  the  1.75  km/s  Impact  of  an  8 
diameter  projectile  against  a  comp-B  target  protected  by  a  4  am  thick  cover 
plate.  In  general,  the  flew  fields  were  quite  similar  to  these  observed  in  the 
bare  charge  cases.  Only  the  case  of  the  6  mm  diameter  projectile  with  the  1.5 
am  thick  cover  plate  showed  the  small  diaaeter  effect  observed  in  the  b.cre 
charge  problems.  In  this  case  we  observed  a  quick  sho^k  to  detonation 
transition  followed  by  quenching.  This  is  shown  in  Figure  13*  However,  we  did 
not  consider  quite  as  small  projectiles  in  the  covered  charge  problem.  With  the 
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Figure  9.  Sequence  of  Mass  Fraction  Contour  Plots  for  the  Supercritical 
Impact  of  an  Annular  Cross-Section  Steel  Projectile  at  0.9  km/ 
against  Bare  Composition-B 
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Figure  9.  (continued) 
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Figure  10.  Sequence  of  Mass  Fraction  Contour  Plots  for  the  Supercritical 

Impact  of  a  "Rectangular'’  Cross-Section  Steel  Projectile  at  1.25  km/s 
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Figure  10.  (continued) 
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Figure  11.  Sequence  of  Mass  Fraction  Contour  Plots  for  the  Subcritical 
Impact  of  a  "Rectangular"  Cross-Section  Steel  Projectile  at 
1.1  ka/s  against  Bare  Corapo3ition-B 


Figure  11.  (continued) 


TIME-  1.05  CYCLE-  210 


r 
10. 

INTERVAL- 


r 

15.  20. 

1 .0000E-01 


25. 


1.25  CYCLE-  250 


Figure  12. 


Sequence  of  Mass  Fraction  Contour  Plots  for  the  Supercritical 
Impact  of  an  8  nun  Diameter  Steel  Projectile  at  1.75  km/s  against 
a  Composition- B  Target  Protected  by  a  4  mm  Thick  Steel  Cover  Plate 
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Figure  12.  (continued) 
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Figure  13.  Sequence  of  Mass  Fraction  Contour  Plots  for  the  Subcritical 

Impact  of  a  6  mm  Diameter  Steel  Projectile  at  1.4  km/s  against 
a  Composition-B  Target  Protected  by  a  1.5  mm  Thick  Steel  Cover 
at 
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Figure  13.  (continued) 
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thickest  cover  plate  (h/d=3/4)  the  rarefaction  was  observed  to  overtake  the 
shock  coapletely  within  the  cover  plate  before  propagating  into  the  explosive. 

In  these  cases,  detonation,  when  produced,  develops  at  a  decaying  shock  wave  and 
the  critical  velocity  is  higher  than  sight  be  anticipated. 

Determination  of  Critical  Impact  Velocity  and  Cosparison  with  Experisental  Data. 
A  United  quantity  of  covered  conp-B  experimental  data  for  comparison  is 
available.  This  includes  the  early  results  of  Slade  and  Dewey  as  well  as  some 

■ore  recent  results  obtained  by  Howe10  for  projectile  attack  against  105  ns 
ammunition.  The  Jacobs-Roslund  formula  for  covered  explosive  suggests  that  the 
product  of  critical  velocity  and  square  root  of  projectile  diameter  depends  only 
on  the  h/d  ratio.  Thus,  in  Figure  14  we  have  plotted  our  2DE  results  together 

1/2 

with  the  aforementioned  experimental  data  in  the  V*d  -h/d  plane.  The  2DE 
predictions  appear  Independent  of  projectile  diameter  but  do  not  produce  a 
straight  line  in  this  plane.  They  agree  quite  closely  with  the  Howe  results  at 
the  two  smaller  h/d  ratios  and  not  as  well  with  the  Slade  and  Dewey  experiments. 
This  is  a  curious  result  since  the  Slade  and  Dewey  experiments,  with  relatively 
thin  cover  plates,  almost  certainly  produce  ahock  initiation  and  Howe  has 
interpreted  his  observed  initiations  with  thicker  shell  casings  as  due  to  a 
shear  mechanism.  The  matter  is  further  complicated  by  the  fact  that  the  2DE 
computations  do  agree  with  the  Slade  and  Dewey  bare  charge  results.  A  straight 
line  has  been  fitted  through  the  2DE  and  Howe  results  for  0.2<h/d<0.6.  Both  the 
experimental  and  theoretical  results  lie  above  the  straight  line  for  h/d>0.6  and 
do  not  agree  closely  with  one  another. 

If  the  2DE  results  for  covered  charges  are  correct,  then  the  Slade  and 
Dewey  results  become  suspect  and  How' *3  experimental  initiations  must  be  due  to 
shock.  The  difference  observed  at  the  larger  h/d  values  would  then  indicate 
lower  computational  accuracy,  possibly  due  to  the  Inadequacy  of  the  Forest  Fire 
model  when  the  rarefaction  immediately  follows  the  initiating  shock.  It  appears 
■ore  likely,  however,  that  the  Slade  and  Dewey  data  are  correct.  Then,  the  2DE 
results  must  be  in  error  and  Howe's  initiations  may  properly  be  attributed  to  a 
mechanism  other  than  shock. 


V.  SUMMARY 

Our  computational  study  of  projectile  impact  shock  initiation  of 
composltlon-B  revealed  details  of  the  flow  fields  produced  and  provided 
predictions  of  critical  impact  velocities  for  both  bare  and  covered  explosive 
targets. 

For  bare  charges,  we  observed  two  different  mechanisms  by  which  the 
critical  velocity  is  determined.  For  Impacts  by  projectiles  of  sufficiently 
large  diameter  initiation  occurs  as  the  impact  induced  shock  wave  builds  to 
detonation  by  reinforcement  due  to  burning  behind  the  shock.  For  smaller 
diameter i  high  velocity  projectiles,  we  saw  that  detonation  or  near  detonation 
breaks  out  immediately  on  impact,  but  may  be  quenched  by  the  enuslng 
rarefactions.  We  found  that  2DE  predicted  the  critical  velocity  accurately.  We 

also  checked  /p2dt  values  along  the  initiation  threshold  and  found  them  to  be 
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relatively  constant.  We  compared  the  shock  to  detonation  transition  paths  to 
the  Pop-plot  for  comp-B  and  found  them  to  agree  in  the  case  of  a  planar  shock 
buildup  but  not  in  the  case  of  projectile  impact,  for  which  multiple  paths  to 
detonation  were  observed.  We  also  simulated  the  special  projectile  geometries 
considered  by  Moulard,  and  found  that  2DE  provided  a  qualitative  explanation  of 
his  observations. 

In  the  case  of  covered  projectiles  we  found  flow  fields  similar  to  the  bare 
charge  case.  The  thickest  cover  plates  allowed  the  rarefaction  to  overtake  the 
shock  before  they  entered  the  explosive  and  significantly  raised  the  critical 
velocity.  The  predicted  initiation  thresholds  agree  with  Howe's  results  but  not 
with  Slade  and  Dewey's. 


v  V.  • 
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3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or 

other  area  of  interest  for  which  the  report  will  be  used.) _ 
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